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CHAOS OF A MARKOV OPERATOR AND THE FOURTH 
MOMENT CONDITION 

By M. Ledoux 

Universite de Toulouse and Institut Universitaire de France 

We analyze from the viewpoint of an abstract Markov operator 
recent results by Nualart and Peccati, and Nourdin and Peccati, on 
the fourth moment as a condition on a Wiener chaos to have a dis- 
tribution close to Gaussian. In particular, we are led to introduce 
a notion of chaos associated to a Markov operator through its iter- 
ated gradients and present conditions on the (pure) point spectrum 
for a sequence of chaos eigenfunctions to converge to a Gaussian 
distribution. Convergence to gamma distributions may be examined 
similarly. 

1. Introduction. In a striking contribution [20], Nualart and Peccati dis- 
covered a few years ago that the fourth moment of homogeneous polynomial 
chaos on Wiener space characterizes convergence toward the Gaussian dis- 
tribution. Specifically, and in a simplified (finite dimensional) setting, let 
F : R N -»• R, 1 < k < N, be defined by 

N 

(1) F = F(x)= ^ a ii,-,ik X h'" x ik' x=(xi,...,x N ) & R N , 
tl,...,tfc=l 

where aj lr .. %, are real numbers vanishing on diagonals and symmetric in the 
indices. Assume by homogeneity that J RJV F 2 d'fN = 1 where 

d lN {x) = {2^)- N ' 2 e -\ x \ 2 l 2 dx 



is the standard Gaussian measure on W . Such a function F will be called 
homogeneous of degree k. Let now F n on R Nn , n E N, N n — > oo, be a sequence 
of such homogeneous polynomials of fixed degree k. The main theorem of 
Nualart and Peccati [20] expresses that the sequence of distributions of the 
F n 's converges toward the standard Gaussian distribution 71 on the real line 
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if and only if 

(2) / F*d lNn ^3 

(3 being the fourth moment of the standard normal). The result actually 
holds for homogeneous chaos on the infinite dimensional Wiener space, and 
the equivalence is further described in terms of convergence of contractions. 
The proof of [20] relies on multiplication formulas for homogeneous chaos 
and the use of stochastic calculus. 

Since [20] was published, numerous improvements and developments on 

this theme have been considered; cf., for example, [13-15, 17, 19, 23], 

An introduction to some of these developments (with emphasis on multi- 
plication formulas) is the recent monograph [22] by Peccati and Taqqu. In 
particular, the work by Nualart and Ortiz-Latorre [19] introduces a techno- 
logical breakthrough with a new proof only based on Malliavin calculus and 
the use of integration by parts on Wiener space. In this work, the conver- 
gence of (-F n ) n6 N to a Gaussian distribution [and thus also (2)] is also shown 
to be equivalent to the fact that 

(3) Var 7jVn (|VF n | 2 )^0, 

where Var 7JVn is the variance with respect to jN n ■ Based upon this observa- 
tion, recent work by Nourdin and Peccati [13, 14] develops the tool of the 
so-called Stein method (cf., e.g., [5, 6, 26, 27]) in order to quantify the con- 
vergence toward the Gaussian distribution. Relying also on multiplication 
formulas and the use of integration by parts on Wiener space, one key step 
in the investigation [13] is expressed by the following inequality: for a given 
homogeneous function F of degree k on normalized in L 2 (7tv), 

(4) Var 7jV (|VF| 2 )<C fc ^F 4 d 7 iv-3), 

where > only depends on k. In particular, the proximity of f RN F 4 djN 
to 3 controls the variance of |V-F| 2 . Now, Stein's method for homogeneous 
chaos on Wiener space as developed in [13] expresses that 

(5) ^T^CVar^lVFl 2 ) 1 / 2 , 

where d(y, 71) stands for some appropriate distance between the law v of 
F and 71, so that |V-F| 2 being close to a constant forces the distribution 
of F to be close to a Gaussian distribution. The conjunction of (4) and (5) 
thus describes how the fourth moment condition controls convergence to a 
Gaussian. 

The primary motivation of this work is to understand what structure of 
a functional F allows for the preceding results, in particular thus the con- 
trol by the fourth moment of the distance to the Gaussian distribution. In 
the process of this investigation, we will revisit the preceding results and 
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conclusions in the setting of a symmetric Markov operator, including, as a 
particular example, the Ornstein-Uhlenbeck operator L = A — x ■ V, corre- 
sponding to the Wiener space setting. In order to achieve this goal, observe 
that the homogeneous polynomial F of (1) is an eigenfunction with eigen- 
value k of the Ornstein-Uhlenbeck operator, that is, —LF = kF. We shall 
therefore try to understand what is necessary for an eigenfunction F of a 
Markov operator in order to satisfy an inequality such as (4). This inves- 
tigation leads us to define a notion of chaos eigenfunction with respect to 
such a Markov operator with pure point spectrum consisting of a count- 
able sequence of eigenvalues, the homogeneous polynomial F of (1) being 
one example with respect to the Ornstein-Uhlenbeck operator. The main 
achievement of this work is then the formulation of an explicit condition on 
the sequence of eigenvalues under which a chaos eigenfunction satisfies an 
inequality such as (4). 

The basic data will thus be a Markov operator L on some state space 
(E,J-) with invariant and reversible probability measure \i and symmetric 
bilinear carre du champ operator 

T(f,g) = ±[L(fg)-fLg-gLf], 

acting on functions /, g in a suitable domain A. For simplicity, we often write 
r(/) = r(/, /) which is always nonnegative. By invariance and symmetry of 
jj, with respect to L, the definition of the carre du champ operator T yields 
the integration by parts formula 

/ f(-Lg)dfi= I g(-Lf)dfi= [ T(f,g)dfi. 
Je Je Je 

In particular J E Lf dfi = since LI = by the Markov property. The oper- 
ator L is said, in addition, to be a diffusion operator if, for every smooth 
function ip : M. — >• R, and every / € A, 

M/) = y/(/)L/ + ^'(/)r(/). 

Alternatively, T is a derivation in the sense that T(ip(f),g) = ip' (f)T(f, g). 

We refer to the lecture notes [1], Chapter 2, by Bakry for an introduc- 
tion to this abstract framework of Markov and carre du champ operators 
and a discussion of some of the examples emphasized below. Additional 
general references include [7] for further probabilistic interpretations and 
[4, 8] for constructions in terms of Dirichlet forms; see also [12] and the 
forthcoming [3]. One prototype example of a Markov diffusion operator is 
the Ornstein-Uhlenbeck operator acting on say the algebra A of polyno- 
mial functions / on E = W N as L/(x) = A/(x) — x ■ V/(x), with invari- 
ant and reversible probability measure the Gaussian distribution /i = 7^ 
and carre du champ = |V/| 2 . One could consider its infinite dimen- 
sional extension on Wiener space (cf. [4] and [18], Chapter 1), but for 
simplicity in the exposition we stick here on the finite dimensional case 
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as a reference example. The preceding general setting also includes dis- 
crete examples, such as the two-point space and its products. Namely, on 
E = {-1, +1}", let Lf = \ Eti Dif where Dif(x) = f{n{x)) - f(x), x = 
(xi, . . . ,Xi, . . . ,xn), Ti(x) = (x±, . . . , — Xi, . . . ,xn). L is invariant and sym- 
metric with respect to the uniform measure n on {— 1,+1} with carre du 
champ r(/) = | J2iLi(Dif) 2 , but is not a diffusion operator. 

These two examples actually entail a crucial chaos structure in the sense 
that the generators L may be diagonalized in a sequence of orthogonal poly- 
nomials (Hermite polynomials in the Gaussian case, Walsh polynomials in 
the cube example); see, for example, [1], Chapter 1, [18], Chapter 1, [10], 
Chapter 2, [22], Chapter 5. More precisely, setting for k = (fei, . . . , kj^) G N N , 
x = (xi, . . . ,x N ) G R N , H k {x) = h kl (xi) ■ ■ -h kN (x N ), with (h k ) keN the se- 
quence of orthonormal Hermite polynomials on the real line, any function 
/ : W N — > R in L 2 (7at) may be written as 

fceN|fc|=fc 

where (-,-) is the scalar product in L 2 (7jv) and where the second sum runs 

over all k £ with \k\ = k\ H h k^ = k. An element H = H k with |fc| = k 

is an eigenfunction of the Ornstein-Uhlenbeck operator with —LH = kH 
and the spectrum of the operator — L thus consists of the sequence of the 
nonnegative integers. For fixed k € N, linear combinations 

(6) F =Y1 a k H k 

\k\=k 

define generic eigenfunctions (chaos) of — L with eigenvalue k, the homoge- 
neous function F of (1) being one example. 
Similarly, if / : {-1, +1}^ R, 

JV 

/=EE (f,w A )w A , 

k=0\A\=k 

where the second sum runs over all subsets A of {1, . . . , N} with k elements 
and 

W A (x) = Y[xi, x = ( Xl ,...,x N )e{-l,+l} N ,Ac{l,...,N}, 

i&A 

are the so-called Walsh polynomials. For the discrete operator L/ = 
h^2i=\Dif , —LW A = kW A if \A\ = k. The spectrum of — L is thus equal 
to N, and linear combinations 

(7) F= *aWa 

\A\=k 

describe the family of eigenfunctions (chaos) of — L with eigenvalue k. 
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A further example is Poisson space. In dimension one, let \x be the Poisson 
law on N with parameter 9 > 0. For a function / :N — > R with finite sup- 
port say, let Df(j) = f(j) - f(j - 1) for every j € N [/(-l) = 0]. The Pois- 
son operator may then be defined as L/(j) = 9Df(j + 1) — jDf(j), j G N. 
It is not a diffusion. The associated carre du champ operator is given by 
2r(/)(j) = 6Df(j + l) 2 +jDf(j) 2 , j G N. The operator — L has a spectrum 
given by the sequence of the integers and is diagonalized along the Charlier 
orthogonal polynomials. Multi-dimensional Poisson models are similar. 

Laplacians L = A on (compact) Riemannian manifolds, and acting on 
families of smooth functions, also enter this framework. These Laplacians 
are diffusion operators and, in the compact case, have again a spectrum 
consisting of a countable sequence of eigenvalues; cf., for example, [9]. 

This work will analyze properties of eigenfunctions of such Markov opera- 
tors L, that is, functions F : E — > R (in the domain of L) such that —LF = XF 
for some A > 0. (We emphasize that F and A are thus rather eigenfunction 
and eigenvalue of — L which is nonnegative.) The ultimate goal of this work 
is to find conditions on such an eigenfunction F of a diffusion operator L in 
order that the analog of (4) holds, and that the fourth moment condition 
then ensures the proximity with the Gaussian distribution. We outline here 
the various steps of the investigation. The first step will be to show (fol- 
lowing [13] in the Ornstein-Uhlenbeck setting) that Stein's method applied 
to an eigenfunction F indicates that it has a Gaussian distribution if (and 
only if) its carre du champ T(F) is constant; see Proposition 1 below. More 
precisely, in accordance with (5), for suitable families of functions tp : R — > R, 
and whenever J E F 2 dfj, = 1 , 



(8) 



ip(F)dfi— / cpdj 



1/2 



where Var^ is the variance with respect to \i. 

On the basis of this result, the fourth moment condition appears quite 
naturally by the integration by parts formula since (assuming the necessary 
domain and integrability conditions) 

A j F 4 dfM= j F 3 (-LF)dfi = 3 j F 2 T(F)dfi. 

Moreover, f E T(F)d/j,= L F(-LF) dfi = A f E F 2 dfj,, so that, still assuming 
by homogeneity that f E F 2 dfi = l, 

(9) AQ^F 4 d/i-l) =J F 2 (T(F) - A) dp. 

This identity is the first indication that the proximity of J E F 4 d/jL with 3 
actually amounts to the proximity of T(F) with its constant mean value A. 

The next step in the investigation, the main result of this note, describes 
a chaos structure of an eigenfunction F of a Markov operator L (not nec- 
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essarily diffusive) with spectrum consisting in a sequence S = {0 = Ao < 
Ai < A2 < • • •} of eigenvalues in order that whenever F is such a chaos with 
eigenvalue A& normalized in L 2 (//), 

(10) Var M (r(F)) < C k [ F 2 (T(F) - X k ) dn 

Je 

for some finite constant C k only depending on 5. The relations (8), (9) 
and (10) together therefore describe how the fourth moment condition 
f E F 4 d[j, ~ 3 ensures that T(F) is close to constant and thus that the dis- 
tribution of F is close to Gaussian. This family of inequalities may then 
be used to describe convergence to a Gaussian distribution of a sequence 
of such chaos eigenfunctions. The abstract chaos structure underlying these 
results is defined by means of the iterated gradients of the Markov operator 
L and is shown to easily cover the examples of Wiener, Walsh or Poisson 
chaos. For example, the chaos structure of the homogeneous polynomial F 
of (6) actually amounts to the fact that V fc+1 i ? = 0. The proof of (10) will 
proceed by a standard and direct algebraic T-calculus on eigenfunctions in- 
volving the iterated gradients of the operator L and avoiding any type of 
multiplication formulas for chaos. 

Turning to the content of this note, Section 2 briefly presents Stein's 
method applied to an eigenfunction of a Markov diffusion operator. The 
next section discusses the iterated gradients and the associated T-calculus on 
eigenfunctions, of fundamental use in the investigation. Section 4 introduces 
the notion of chaos of a Markov operator with pure point spectrum and 
presents the aforementioned main result (10), proved in Section 6. The last 
section briefly describes analogous conclusions for convergence to gamma 
distributions covering recent results of [14]. 

It should be carefully emphasized that the present exposition develops 
more the algebraic and spectral descriptions of the problem under investiga- 
tion [and concentrates on a proof of (10)] rather than the analytic issues on 
domains and classes of functions involved in the analysis. In particular, we 
work with families of functions in the domain of the Markov operator and its 
carre du champ and with eigenfunctions assumed to satisfy all the necessary 
domain and integrability conditions required to develop integration by parts 
and the associated T-calculus. These properties are classically and easily sat- 
isfied for the main examples in mind, the Gaussian case, the discrete cube 
or the setting of the Laplace operator on a compact Riemannian manifold. 
Note, however, that the extension from the finite dimensional Gaussian set- 
ting to the infinite dimensional one requires basic analysis on Wiener space 
as presented, for example, in the first chapter of [18] (see also [22]) in order 
to fully justify the domain issues and the various conclusions. These aspects, 
carefully developed in the aforementioned references, are not discussed here. 
Further conditions ensuring the validity of the results presented here might 
be developed in broader contexts. 
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2. Stein's method for eigenfunctions. We start our investigation with 
a brief exposition of Stein's lemma applied to eigenfunctions of a diffusion 
operator. We refer to [5, 6, 26, 27] and the references therein for general 
introductions on Stein's method. The results below are mere adaptations 
of the investigation [13] by Nourdin and Peccati in Wiener space to which 
we refer for further details. Throughout this section, L is thus a diffusion 
operator with invariant and reversible measure \i and carre du champ T as 
described in the Introduction. All the necessary domain and integrability 
conditions on the eigenfunctions under investigation are implicitly assumed, 
and are satisfied for the main Ornstein-Uhlenbeck example; cf. [13]. 

We first illustrate, at a qualitative level, Stein's method in this abstract 
context. Given a measurable map F:E—>M., say that L commutes to F 
if there exists a Markov operator C on the real line such that for every 
tp : R — > R (in the domain of C and such that <p o F is in the domain of L) 

L(<poF) = (&p)(F). 

In this case, the image measure \ip of fi by F is the invariant measure of C 
One model factorization operator C on R is the Ornstein-Uhlenbeck op- 
erator Ctp = tp" — xtp' with invariant measure the standard Gaussian dis- 
tribution d'ji(x) = e~ x I 2 tt|=- Let then F be an eigenfunction of — L with 
eigenvalue A > 0. The observation here, at the root of Stein's argument, is 
that whenever F = F(F) is (^-almost everywhere) constant, then L com- 
mutes to F through the Ornstein-Uhlenbeck operator C, and thus the dis- 
tribution fip of F is Gaussian. Namely, note first that by integration by 
parts, f E Fdfi = J E F(—LF) dfi = A f E F 2 d/u, so that if F is constant and F 
is normalized in L 2 (/i), then F = A. Then, for ip : R — > R smooth enough, by 
the chain rule formula for the diffusion operator L, 

L(ip o F) = ip'(F)LF + (p"(F)F = -XFip'(F) + <p"(F)F. 
Hence, if F = A, 

L(ipoF) = \(£<p)(F) 

so that L commutes to F, and thus /xp is the invariant measure of the 
Ornstein-Uhlenbeck operator C characterized as the Gaussian distribution 71 
For an eigenfunction F, F = F(F) constant thus forces the distribution 
of F to be Gaussian. Now, as such, this observation is not of much use 
and to describe convergence to normal as for sequences of homogeneous 
polynomials in the Introduction, it should be suitably quantified in the form 
of inequality (8) in order to express that the proximity of F with a constant 
value forces the distribution of F to be close to Gaussian. This is the content 
of the classical Stein lemma as described in the next statement. 

Proposition 1. Let F be an eigenfunction of — L with eigenvalue A > 
and set F = F(F). Denote by \xp the distribution of F. Given (p :R —> R inte- 
grable with respect to [ip and 71, let tp be a smooth solution of the associated 
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Stein equation ip — J^cpdji = tp' — xtp. Then, 



ill) 



[ ydii F - f <pd 7l <^f([ (T-X) 2 dfi) 



1/2 



where = H^'lloo- ^ n particular, if f E F 2 d[i = 1, 



cpdfj,j 



<^Var M (T)V». 



Proof. Since hf is the distribution of F under fx, and by the Stein 
equation, 



(pdfj,} 



<p{F)dfi- [ ^d 7 i= I bP\F)-FiP(F)]d(i. 

JR JE 



Now —LF = XF so that 

ip'(F) - Ftp(F) = tp'(F) + X^LFipiF) 

and hence, after integration by parts with respect to the operator L and the 
use of the diffusion property, 

I ipdfip- [ ipd ll = I ^'(F)[l-X- l T)dfi. 

JR JR JE 



Together with the Cauchy-Schwarz inequality, 

1/2 



ipdfii 



ipdji 



< 



Jj\Ffd^\ ^j[l-X- l Tfd^ 



1/2 



which amounts to (11). If f E F 2 dfi = l, then J E T dfi = f E F(—LF)dfi = X 
and thus f E {T — A) 2 d\i = Var A1 (r). The proof of Proposition 1 is complete. 

□ 



Proposition 1 is thus investigated in [13] for Wiener chaos. As is discussed 
there (Lemma 1.2 and Theorem 3.1), the constant in (11) of Proposition 1 
can be uniformly bounded inside specific classes of functions. For instance, 
Cp < 2 when ip is the characteristic function of a Borel set (corresponding 
to the total variation distance) and C v < 1 when ip is the characteristic 
function of a half- line (corresponding to the Kolmogorov distance). 

For the further purposes, observe, as is classical (cf. [26, 27]), that Stein's 
strategy may be developed similarly for the Laguerre operator on the positive 
half-line C p ip = xip" + (p — x)ip' ' , p > 0, with invariant measure the gamma 
distribution dg p (x) = r(p)~ 1 x p ~ 1 e~ x dx. Let F be an eigenfunction of — L 
with eigenvalue A > and V = T{F). As above, for every ip : IR — > IR smooth 
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enough, setting G = F + p, 

L(<p oG) = <p'(G)LF + ip"{G)T 
= -\F<p , (G)+<p"(G)T 



xf(p-G)^(G) + ^"(G)\ 



In this case, if T = XG, 

%oG) = A(£ pV )(G) 

so that hg is the invariant measure of C p characterized as the gamma dis- 
tribution g p . 

For this example of the Laguerre operator, the criterion for an eigenfunc- 
tion F to have a gamma distribution is thus that T = X(F + p). On the basis 
of this qualitative description of Stein's method for the Laguerre operator, 
the next statement illustrates the analog of Proposition 1 for this model. 

Proposition 2. Let F be an eigenfunction of — L with eigenvalue X > 0, 
and set T = T(F). Let p > and denote by hf+p the distribution of F + p. 
Given (p:~R— >M, integrable with respect to ^f+p and g p , let tp be a smooth 
solution of the associated Stein equation (p — J" R ip dg p = xtp 1 + (p — x)ip. Then, 

1/2 

(12) 



/ tpdfiF+p - [ ipdg p <%f/ (r - X(F + p)) 2 d[i) , 

JR JR A \JE J 



where C a 



. In particular, if J E F 2 d[i = p, 



ipdg % 



^Var^r-AF) 1 / 2 . 



Proof. Set again G = F + p. Start as in the proof of Proposition 1, 
namely 

(pd/j, G - / fdg p = / (p(G)dfi- / <pdg p 

JR JE JR 

[G^(G) + (p-G)i;(G)]d f i. 



Since —LF = XF, and thus LG = X(p - G), 

Gip'(G) + {p- G)ip{G) = Gip'{G) + A _1 LGy>(G). 

After integration by parts with respect to the operator L and the use of the 
diffusion property, 

f vdfiG- f <pdg p = [ ip'(G)[G-X- x T]dLi. 

JR JR JE 

The conclusion follows similarly from the Cauchy-Schwarz inequality. □ 
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Proposition 2 is similarly investigated in [13] in the context of Stein's 
method on Wiener space. Again the the constant in (12) may be bounded 
only in terms of p inside specific classes of functions; cf. [13], Lemma 1.3 
and Theorem 3.11. Analogs of Stein's lemma in the context of the preceding 
statements have been investigated on discrete Poisson or Bernoulli spaces 
in [16, 21, 24]. In those examples, the control of the variance of T is not 
enough to ensure proximity to a Gaussian distribution and has to be sup- 
plemented by various additional conditions. 

3. Iterated gradients. This section presents the family of the iterated 
gradients of a Markov operator and the basic (algebraic) T-calculus on eigen- 
functions at the root of the investigation. Given a symmetric Markov opera- 
tor L as above (not necessarily a diffusion operator), recall following [1, 11], 
the iterated gradients L m , m>2, associated to L defined according to the 
rule defining T = Ti as 

T m (f,g) = i[Lr m _i(/ )5 ) -r m _i(/,L 9 ) -r m _i( 5 ,L/)] 

for functions f,g in a suitable class A. By extension, To(f,g) = fg. For 
simplicity, set T m (f) = T m (f, /). Note that in general T m (f) for m > 2 is not 
necessarily nonnegative. The T2 operator has been introduced first by Bakry 
and Emery [2] to describe curvature properties of Markov operators and to 
provide a simple criterion to ensure spectral gap and functional inequalities; 
cf. [1], Chapter 6, [12] and [3]. This criterion will be used in Proposition 4 
below. The iterated gradients T m have been exploited in [11] toward variance 
and entropy expansions. 

The following elementary lemma will be of constant use throughout this 
note and concentrates on the significant properties of the iterated gradients 
of a given eigenfunction. Recall that we assume the necessary domain and 
integrability conditions to justify the relevant identities. 

Lemma 3. Let F be an eigenfunction of — L with eigenvalue A. Set T m = 
T m (F), m> 1. Then, for every m>\, 

(13) r m = iLr m _! + Ar m _x = (±l + a id) m_1 r. 

Furthermore, for every m,n> 1, 

(14) / T n T m dfi= / r n _ir m+ id/i. 
J E Je 

In particular, by selecting n = 1 , for every m > 1 , 



(15) 



f TT m dfi= [ F 2 T m+1 dfi. 
Je Je 
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Proof. Equality (13) is an immediate consequence of the definition of 
r m and the eigenfunction property 

r m (F) = iLr m _i(F) -r m _i(F,LF) = ±Lr m _!(F) +Ar m _!(F). 

The conclusion follows by iteration. 

Recalling the notation T m = F m (F), multiply the preceding identity by 
T n and integrate with respect to /x to get, by symmetry, 

2 / r n r m dfi= / r m _iLr n (i/i + 2A / r n r OT _id/x. 

Je Je Je 

Changing the role of n and m — 1, by symmetry again, 

2 / r m _ir n+ i dfi= / r m _iLr n d/i + 2A / r m _ir n dp 

Je Je Je 

and the identity (14) follows. The proof of the lemma is complete. □ 

The following statement is a first illustration of the method developed 
next. It expresses a kind of rigidity result under the geometric T2 curvature 
condition mentioned previously. 

Proposition 4. Assume that the operator L is of curvature p > in 
the sense of Bakry-Emery [2] ([1], Chapter 6), that is, I^/) > /°P(/) f or 
every f 6 A. If F is an eigenfunction of — L with eigenvalue p, then T(F) 
is (n-almost everywhere) constant. In case L is a diffusion operator, the 
distribution of F is Gaussian. 

It might be useful to recall ([1], Chapter 6, [3, 12]) that under the curva- 
ture condition of the statement, A > p for every nonzero eigenvalue A of — L. 
In particular, L is ergodic in the sense that if T(/) = 0, then / is constant 
(//-almost everywhere). It is also worthwhile mentioning that for the model 
space consisting of the Ornstein-Uhlenbeck diffusion operator L = A — x ■ V 
with invariant measure 77V , p = l and the eigenfunctions with eigenvalue 1 
are the linear functions 

N 

F(x) = ^aiXi, x = (x 1 ,...,x N ) £R N , 
i=i 

whose distributions are of course Gaussian. Since Gaussian Wiener chaos 
of order larger than or equal to 2 do not contain any nonzero Gaussian 
variable [10] and [20], Proposition 4 thus expresses a kind of rigidity property 
in the sense that if F is a nonzero eigenfunction of the Ornstein-Uhlenbeck 
operator L with eigenvalue A, then F is Gaussian if and only if T{F) is 
constant, and if and only if A = p = 1. 

The proof of Proposition 4 is rather straigthforward. Write as before 
T m (F) = T m , m > 1. By Lemma 3 [formula (13)], T2 = ^Lr + pT. There- 
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fore, under the curvature condition ^(J) > pT(f), LP > 0. But then 

0< / TLTdfj, = - / T(T)dfi<0, 
Je Je 

so that r = T(F) is (//-almost everywhere) constant. The final assertion of 
the statement then follows from Stein's lemma (Proposition 1). 

4. Chaos of a Markov operator. This section is devoted to the main 
conclusions of this work. We are thus given, on a state space E, a Markov 
operator L with symmetric and invariant probability measure /i and carre du 
champ r (acting on a suitable algebra of functions A). Assume in addition 
that L has a pure point spectrum consisting of a countable sequence of 
eigenvalues S = {0 = Ao < Ai < A2 < • • •} (more precisely, S is the spectrum 
of — L) (cf. [3, 25, 28]). Since Ai > 0, L is ergodic [in the sense that if T(/) = 0, 
then / is constant]. 

Given the spectrum 5 = {0 = Ao < Ai < A2 < • • •}, define for every k G N 
the polynomial of degree k in the real variable X, 

fc-i k 

Q k (x) = \{(x-x l )=Y J ^Qfm t 

i=0 i=l 

(Qo = 1). Define then the bilinear form (acting on Ax A) 

k 

i=i % - 

The following main definition introduces the notion of chaos associated to L 
and its spectrum S. 

Definition 5. An eigenfunction F of — L with eigenvalue \ k (— LF = 
XkF) is said to be a chaos of degree k > 1 relative to S if Q k+ x(F)(F) = (/i- 
almost everywhere). We call F a chaos eigenfunction (with eigenvalue Afc). 

Motivation for the preceding definition is provided by the Ornstein-Uhlen- 
beck operator with spectrum S = N. Namely, it is easily shown in this case 
(see [11], Section 2) that Q k (T)(F) = \V k F\ 2 . Any eigenfunction F as in (6) 
is such that V k F is constant and V k+1 F = leading thus to Definition 5. 
In the infinite dimensional setting of an abstract Wiener space (E,H,n) 
with separable Hilbert space H, referring to [18], Chapter 1, for notation 
and terminology, the Ornstein-Uhlenbeck operator L has domain B 2 ' 2 and 
Qfc(r)(F) = ||D fc F||^® fe for any F G D fc ' 2 where D is the derivative opera- 
tor (use as in the finite dimensional case the commutation [L,Z?] = D and 
the chain rule formula [18], Proposition 1.4.5). Now, if J k F denotes the 
projection of F (in D fe ' 2 ) on the fcth Wiener chaos, LJ^F = —kJ k F and 
D k (J k F) = J D k F = E(D k F) so that J k F thus defines a /c-chaos in the 
sense of Definition 5. For example, in case H = L 2 (T, B, v) where v is a o- 
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finite atomless measure on a measurable space (T, £>), the elements J k F may 
be represented as multiple stochastic integrals 

Wk)= I ■■■ [ fk(h,...,t k )w(dt 1 )---w(dt k ) 

JT JT 

of symmetric functions f k on L 2 (T fc ) with respect to the white noise W and 
D k I k {f k ) = {f k (h,. . .,t k );h,. . . ,t k £ T}. 

The discrete operator L/ = iXa^i-^*/ 011 the CUDe { — Ij+I}^ an d the 
Poisson operator are further instances entering this definition with again 
5 = N (see [11], Section 2). On the cube {— 1,+1} , for example, 

Qk(T)(F) = ^^ D ^--- D ^ 2 ' 

where the sum is over distinct i\, . . . ,i k £ {1, . . . , N} and thus any F of the 
form (7) is a A;-chaos (k < N). 

There are of course examples of eigenfunctions which are not chaos. For in- 
stance, the Laguerre operator on the positive half-line C p ip = xip" + (p— x)ip' , 
p > 0, has spectrum equal to N (with eigenvectors the Laguerre orthogonal 
polynomials with respect to the gamma distribution g p ), but the eigenfunc- 
tion F = x — p with eigenvalue 1 is not a 1-chaos as Q2{T)(F) = —\F. 

According to the preceding examples, another possible definition of k- 
chaos would have been that Q k (T)(F) is constant. (If F is normalized in 
L 2 (n), then ([11], page 443), 

I Q k (T){F)dp = [ FQ k (-L)Fd^ = Q k (X k ), 

JE J E 

hence Q k (T)(F) = Q k {\ k ) ) Now, it is easily checked [using (13) of Lemma 3] 
that if F is an eigenfunction of — L with eigenvalue then LQ k (T)(F) = 
2Q k+ i(T)(F). In particular therefore, if Q k (T)(F) is constant, then 
Q k+1 {T)(F) = 0. Conversely, if Q k+1 (T)(F) = 0, by ergodicity, Q k (T)(F) is 
constant. It will turn out more simple in the proofs of the main results to 
use the first definition of chaos [as Q k+ i(T)(F) = 0]. 

The following statements are the main results of this work. Recall the 
polynomials Q k (X) and set, for k > 1, X £ R, 

R k+l (x) = ^\Q k+l (x) - q£>i(o)*] = E lQkU°) xi ~ 2 

i=2 

and 

T k+ \{X) = R k+ i(X + Afc) — R k+ i(X k ). 

Thus, for example, Q 2 (X) = X 2 - X 1 X, R 2 = 1 and T 2 = 0, Q 3 {X) = X 3 - 
(Ai + \ 2 )X 2 + A1A2X, R 3 (X) = X - (Ai + A 2 ) and T 3 (X) = X. Set further- 
more 



Tr k = X 1 ---X k , k>l (vr = l). 
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The following theorem puts forward the fundamental identity at the root 
of this work. 

Theorem 6. In the preceding setting, let F be a k-chaos eigenfunction 
with eigenvalue X k , k > 1. Set T = T(F). Then 



(16) 7T k 



s _! Jj 2 dfi = 7T k J^ F 2 T dfi + (-l) k TT k+1 H^j r dfx. 



Corollary 7. In i/ie preceding setting, let F be a k-chaos eigenfunction 
with eigenvalue X k , k > 1. Sei T = T(F). If 

(17) (_i)fe Tfe+ /-M <0 for every n € N, 
then 

(18) / r 2 d^< A fc f F 2 Tdfi. 

Je Je 

In particular, if F is normalized inL 2 ([i), then f E Tdfi = f E F(— LF) dpi = 
Afc and thus 



Var M (r)<A fc ^F 2 r, 



(19) Var M (r)<A fc / F 2 T dfx - \ k 



Under the additional diffusion hypothesis on L, according to (9), inequal- 
ity (19) of Corollary 7 may be expressed equivalently as 



(20) 



Var^r^AlQyV^u-l). 



In particular, if f E F 4 dfi = 3, then V = V(F) is constant and by Stein's 
lemma (Proposition 1), the distribution of F is Gaussian. 

The next statement describes a fundamental instance for which the spec- 
tral condition (17) in Corollary 7 is fulfilled. 

Theorem 8. The spectral condition (17) in Corollary 7, 
{-l) k T k+1 (-YJ < for every n G N 
is satisfied when S = (A n ) ngN = N. 

As a consequence of this result, the conclusions of Corollary 7 apply to 
the examples of the Ornstein-Uhlenbeck, Bernoulli and Poisson operators. 
As such, some of the main conclusions of [13] are covered by the preceding 
general statement, and in particular the initial result of [20], namely that 
if (F n ) n£N is a sequence of homogeneous Gaussian chaos, normalized in 
L 2 (7Ar n ), N n — > oo, then (F n ) ngN converges to a Gaussian distribution as 
soon as f E F^dfj,—)- 3. 
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For discrete models as the cube or the Poisson space, the picture is less 
satisfactory. For instance on the cube E = {— 1,+1} ™, N n — > oo, if F n = 
X^|A|=fc aWA) n G N, is a sequence of Walsh chaos of degree k normalized 
in L 2 (/i) for the uniform measure fi, and if f E F 2 r(F n ) d/j, — > k, then as an 
application of Corollary 7, T(F n ) — > k in L 2 (fj,). Now T(F) being constant in 
this case is not always discriminative [as shown by the example of F(x) = 
x\ ■ • ■ Xk] and further conditions have to be imposed on the sequence (F n ) neN 
to ensure convergence toward a Gaussian distribution. This analysis has been 
recently achieved in [16]. Similar additional conditions have been studied on 
Poisson spaces in [21, 24]. The input of Corollary 7 on convergence of chaos 
in these discrete examples is that it reduces the convergence T(F n ) — > Xk in 
L 2 (//) by the weaker condition f E F 2 T(F n ) d\x — > Xk- 

5. Chaos of order 1 and 2. Before turning to the general proofs of The- 
orem 6 and Corollary 7, and to get a better feeling about these statements, 
we discuss in this section the particular values k = 1 and k = 2. Recall that 
we write for simplicity T m = T m (F), m > 1, for an eigenfunction F. 

When k = 1, that is, Q 2 (T) = T 2 - XiT = 0, multiplying this identity by F 2 
and integrating with respect to fi, it follows thanks to Lemma 3 [formula (15)] 
that 

/ T 2 dfi = X 1 [ F 2 Tdfi. 
Je Je 

Now here R 2 = 1, and thus T 2 = 0, so that both the fundamental identity (16) 
and the spectral condition (17) are automatically satisfied. 

When k = 2, start from Q 3 (T) = T 3 - (Ai + A 2 )r 2 - AiA 2 r = 0. Multiplying 
by F 2 and integrating, it follows similarly thanks to Lemma 3 [formula (15)] 
that 

/ rr 2 dfi - (Ax + A 2 ) I T 2 d^ + Ai A 2 [ F 2 T dfi = 0. 
Je Je Je 

By (13) of Lemma 3, T 2 = ^Lr + A 2 r so that 

- I rLr d/j, — Ai f T 2 dfi + X l X 2 [ F 2 Tdfi = 0. 
2 Je Je Je 

Here R 3 (X) = X — (Ai + A 2 ) and T 3 (X) = X so that the fundamental iden- 
tity (16) holds, and the spectral condition (17) amounts to X n > for every 
n e N. 

One observation on which we will come back in the next section is that, 
in the case k = 2, only the inequality Qs(T) > is used in order to reach the 
conclusions of Corollary 7. A further observation is that for chaos of order 
1 or 2, the spectral condition (17) is fulfilled for any sequence of eigenvalues 
= Aq < Ai < A 2 < • • • . This is clearly not the case when k > 3. 
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6. Proofs of Theorems 6 and 8. In this section, we establish Theorem 6, 
Corollary 7 and Theorem 8. Let thus F be a fc-chaos with eigenvalue A^. (If 
necessary, we may assume that k > 3 according to the preceding section.) 
Write as usual T m for T m (F), m > 1. 

As in the preceding section for chaos of order 1 or 2, start as a first step 
from the chaos hypothesis Qk+i(P) = 0. Multiply this identity by F 2 and 
integrate with respect to fi. By definition of Qk+i and (15) of Lemma 3, 

0= / F 2 Q k+1 (T)dv 

JE 

(21) = E^ +1 (0) / F 2 T % d^ 

i=l *• JE 

Je i=2 *i Je 

Now, by (13) of Lemma 3, 
fc+i 

Ettq£.i(o) / rr «-i^ 

i=2 l - JE 

fc+i r f \ \ i ~ 2 

= E^l+i(°) / r -L + A fc id rd M 

i=2 ' Je \ / 

=E^ (o) 2C/)^ 2 "/ rL£r ^ 

j=2 ' £=0 V 7 ^ 

= E£ ( i < 2 )^>)>r 2 -^/rL'r< iM . 

£=0i=f+2 V / - JE 

Recalling the definition of the polynomial Rk+i, note that 
fc+i 

e (VH^«< o > A *~ 2 ~'4^*>- 

i=£+2 V 7 

Hence 

fc+i fc— l 

E^fciiW / rr,-i^ = E^£i(^)i / rL ' r ^- 

i=2 ' JE 1=0 ' JE 

Now 

fc-i 

E Tj^fc+it W = R k+1 (X + A fe ) = T k+1 (X) + R k +i(Xk) 
i=o l - 
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so that 

YsWkliW I rr Vi^= / TT k J^\YdpL + R k+1 {X k ) ! T 2 dfi. 

■ 2 * J E J E \ / J E 

The fundamental identity (16) of Theorem 6 then follows from (21) together 
with the fact that 

Ql : 1 i(0) = (-l) fe A 1 ---A ifc = (-l) fe 7r fe 

and 

R k+ i(\k) = (-l) fc+1 Ai • • • A fc _! = (-l^+V-i- 

The proof is complete. 

Corollary 7 is deduced from Theorem 6 through the following classical 
and elementary property, consequence of the point spectrum hypothesis. 

Lemma 9. If P is a polynomial, J E uP(L)udfi > for every u [in the 
L 2 (fi)- domain of P(L)] if (and only if) P(—X n ) > for every n € N. 

Proof. For each n G N, denote by E n the eigenspace associated to the 
eigenvalue A n so that L 2 (/i) = © ngN -E n since S = (A n ) ngN is the spectrum 
of L. Decompose then u in L 2 (//) as u = X^neN n « with u n € E n , n € N, so 
that P(L)u = YlneN P (~^)un and 



f uP{V ) udii = Y J P{-K) [ u 2 n 



dfi 



from which conclusion follows. □ 



As mentioned for chaos of order 2, when k is even, only the inequality 
Qfc + i(r) > is used in order to reach the conclusions of Corollary 7. 

We next turn to the proof of Theorem 8, checking the spectral condition 
(17) {-l) k T k+l (-^f) < 0, n G N, for S = (A n ) n6N = N. Since in this case 
T k+1 (X) = R k+ i{X + k) - {-l) k+1 {k - 1)!, we have to show that 



2 

2 



When ^ = k, the expression on the left-hand side is equal to k\ J2i=i J so 
that the conclusion holds in this case. When ^ ^ k, we need to show that 



> (k - 1)1. 
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Assume first that n > 2k + 1. Then 

n(H=(i-* +1 )n(i 

i=0 v 7 v 7 i=0 v 

A- 



*(i-* +1 )nHM*-* +i ) H - 



Hence 



x -i rfe-i / \ 

) S 5-« -« 



>*!, 



.i=0 

which answers this case. We turn to the case where n < 2k — 1 for which it 
is necessary to check that 

k-l 



i— n \ / 



i=0 

It is enough to assume that n is odd, = 2p — 1, 1 < p < fc. Then 

k-l , x p-1 / x fc-1 / x p , x fc-p / 

n (H=n (i-On f-» ^n -- n - 

i=0 v 7 i=0 v 7 i=p v 7 i=l v 7 i=l v 

Therefore, the inequality to establish amounts to 

n (<-2 n (i-^)<(p-mk-p)i<(k-iy., 

i=l ^ 7 i=l ^ 7 
which is trivially satisfied. The claims thus holds in this case too. Theorem 8 
is therefore established. 

7. Convergence to gamma distributions. In this last section, we briefly 
address the analogs of Theorem 6 and Corollary 7 in the context of con- 
vergence to gamma distributions on the basis of the corresponding Stein 
characterization of Proposition 2. The main conclusion is obtained by a 
simple variation on the fundamental identity (16) of Theorem 6. In particu- 
lar, the analysis covers the recent results of [14] (see also [13]) in the context 
of Wiener chaos. 

The framework is the one of the preceding sections, with a Markov oper- 
ator L with spectrum S = (A n ) ngN and invariant and reversible probability 
measure fj, and carre du champ T. Recall nk = M ■ ■ ■ k > 1, and the poly- 
nomials Rk+i and Tfc+i of Theorem 6. 

The following theorem addresses approximation of a fc-chaos F by a 
gamma distribution via the control of Var^(r — XkF) as emphasized in 
Proposition 2. As announced, the proof is an easy modification on the fun- 
damental identity (16) of Theorem 6. 
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Theorem 10. Let F be a k-chaos with eigenvalue X k , k > 1, such that 
j E F 2 dfi = p > 0. Set T = r(i ? ). Under the spectral condition (17) 
(-l) fc T fc+1 (-4f ) < for every n G N, it holds 

Nm^T -X k F)<X k / F 2 Tdy + A k [ FT dy — pB k — p 2 X 2 ., 
Je Je 

where 

TTfc-l V 2 J ^k-l V 1 

In the diffusion case, 

X k / F 4 dfi = 3 / F 2 Tdy and A fc / F 3 dfi = 2 / Frd/x 
j_e Je Je Je 

so that the conclusion of the theorem reads 

Var,(r-A fc F)<f [ F Uy + ^ [ F 3 dfi —pB k —p 2 X\. 
•j Je 1 Je 

Consider now the example where S = N for which we know from Theo- 
rem 8 that the spectral condition (17) holds. The inequality of Theorem 10 
takes a nicer form when k > 2 is even. Indeed in this case (— l) k X k R k+ i(^-) = 
—2k\ so that 

|var M (r- X k F) < f F 2 T dfi - A f FT dfi + 2pk —p 2 k. 
k J E Je 

In particular in the diffusion case, 

(22) -prYaiJT - X k F)< [ F 4 dfi - 6 [ F 3 dfi + 6p - 3p 2 . 

k 2 Je Je 

This inequality (22) then ensures, through Stein's lemma (Proposition 2), 
that if (-F n ) rie N is a sequence of fc-chaos such that j E F 2 dfi = p for every n 
and 



/ F 4 dfi-6 I F 3 dfi + 6p-3p z -> 0, 
Je Je 



then [F n +p) ng pj converges in distribution to the gamma distribution with 
parameter p, that is the main result of [14]. 

Proof of Theorem 10. Let thus F be a A;-chaos with j E F 2 dfi = p, 
hence f E T dfi = pX k . Set U = T — X k F (so J E U dfi = pX k ). It is immediately 
checked that 

I T 2 dfi = [ U 2 dfi + 2X k [ FTdfi-pX 2 k 
Je Je Je 

= Var M ([/) + 2A fe / FT dfi - p(l - p)X 2 k 
Je 
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and, for every I > 1, 

f TL e T dfi = [ UL e Udfi + 2{-lY\ e k +1 [ FTdn-p(-lYX e + 2 . 
Je Je Je 

Therefore, the fundamental identity (16) of Theorem 6 takes the form, after 
a little algebra, 

(-l) fc J UT k+1 {^jUdu - TTfe-! Var M (C7) 

+ 7r k jT F 2 r dfi + 2(-l) fc A fc i4+i ( ^ J J FT dpi 

-p(-lf\lR k+l (^\ -p 2 X 2 7 r k _ 1 = 0. 
Under the spectral condition (17) (— l) k T k+ i(— 4p) < for every n € N, 
7r fc _iVar M (l7) <ir k J^F 2 Tdv + 2(-l) k \ k R k+1 (^J J FT dfi 

which amounts to the statement of the theorem. The proof is complete. □ 
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